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Abstract 

This paper presents a study on the spectral evolution of gamma-ray burst (GRB) prompt emissions 
observed with the Suzaku Wide-band All-sky Monitor (WAM). By making use of the WAM data archive, 
6 bright GRBs exhibiting 7 well-separated fast-rise-exponential-decay (FRED) shaped light curves are 
presented and the evaluated exponential decay time constants of the energy-resolved light curves from 
these FRED peak light curves are shown to indicate significant spectral evolution. The energy dependence 
of the time constants is well described with a power-law function t{E) oc E 7 , where 7 ~ —(0.34 ±0.12) 
in average, although 5 FRED peaks show consistent value of 7 = —1/2 which is expected in synchrotron 
or inverse-Compton cooling models. In particular, 2 of the GRBs were located with accuracy sufficient 
to evaluate the time-resolved spectra with precise energy response matrices. Their behavior in spectral 
evolution suggests two different origins of emissions. In the case of GRB 081224, the derived 1-s time- 
resolved spectra are well described by a blackbody radiation model with a power-law component. The 
derived behavior of cooling is consistent with that expected from radiative cooling or expansion of the 
emission region. On the other hand, the other 1-s time-resolved spectra from GRB 100707A is well 
described by a Band GRB model as well as with the thermal model. Although relative poor statistics 
prevent us to conclude, the energy dependence in decaying light curve is consistent with that expected in 
the former emission mechanism model. 

Key words: gamma rays: burst — gamma rays: observation — radiation mechanisms: thermal - 
blackbody radiation; non-thermal — inverse Compton scattering 



1. Introduction 

Gamma-ray bursts (GRBs) arc the most energetic 
events known to occur in the universe. Their luminous 
yet variable behavior requires efficient radiation from rel- 
ativistic outflows, although observations have not yet al- 
lowed us to obtain a clear understanding of the radia- 
tion process. Prompt emission spectra have been observed 
mostly in the gamma ray band and have been described 
with a power law, an exponentially cut-off power law, or 
a smoothly connected broken power law (GRBM) (e.g., 
Band et al. 1993; Preece et al. 1998; Preece et al. 2000; 
Barraud et al. 2003; Kaneko et al. 2006; Kaneko et al. 
2008; Abdo et al. 2009; Gonzalez et al. 2009). These 
non-thermal models are consistent with the energy spec- 
tral distributions from optically thin synchrotron or syn- 
chrotron sclf-Compton radiation produced by relativistic 
electrons accelerated through shock fronts in the outflows 
(Tavani 1996; Cohen et al. 1997; Schaefcr et al. 1998; 
Frontera et al. 2000; Wang et al. 2009). 



An energy peak seen in the cut-off power law or the bro- 
ken power law is then attributed to the low energy end of 
the electron energy distribution. However, although non- 
thermal models have succeeded in describing the time- 
averaged spectra of prompt emissions, some observational 
results seem to be inconsistent with the models, such as 
in low-energy spectra exhibiting flatter slopes than those 
allowed in the optically thin synchrotron or synchrotron 
sclf-Compton (SSC) model obtained in some time-resolved 
spectral analyses (Crider et al. 1997; Preece et al. 2002; 
Ghirlanda, Celotti, & Ghisellini 2003; Gonzalez et al. 
2009; Ohno et al. 2009; Onda et al. 2010). These findings 
have motivated researchers to study other emission mech- 
anisms, such as small-pitch-angle synchrotron emission 
(Lloyd & Petrosian 2000; Medvedev 2006) and inverse- 
Compton (IC) scattering of a self-absorbed seed emis- 
sion, although in these cases contradictions arise from the 
point of view of acceleration of electrons (Spitkovsky 2008; 
Giannios & Spitkovsky 2009), inefficiency of energy release 
(Kobayashi et al. 1997; Lazzati et al. 1999; Nysewandcr 
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Table 1. Event list 



WAM trigger 


GRB ID 


Trigger time (UT) 


Position (J2000.0) 


T90 


Reference* 


# 0058 


050924 


08:49:12 


not determined 


8.41 ±0.05 s 


(1) 


# 0408 


060915 


08:25:34 


not determined 


43. 15 ± 0.07 s 


(2) 


# 0415 


060922 


17:21:21 


not determined 


13.88 ±0.26 s 


(3) 


# 0449 


061101 


21:26:28 


not determined 


21.57±0.02 s 


(4) 


# 1129 


081224 


21:17:55 


(14 h 12 m 35?85, ±74 d ll m 59f5) 


11.67±0.20 s 


(5) 


# 1555 


100707A 


00:46:36 


(15 h 03 m 52?41, -30 d 39 m 30?2) 


9.88 ±0.41 s 


(6) 


* (1) Konus GRB #2634; (2) Konus GRB #2790; (3) Konus GRB #2796 / SPI-ACS; (4) Konus GRB #2808; (5) GCNC 8739; (6) 



GCNC 10950/10947 



et al. 2009) or the observed relation between luminosity 
and energy peaks (Ramirez-Ruiz & Lloyd- Ronning 2002). 

Alternatively, an optically thick thermal (blackbody) 
component has been proposed in addition to the pure 
non-thermal (power-law) model (Meszaros & Rees 2000; 
Meszaros et al. 2002; Daigne & Mochkovitch 2002; Rees 
& Meszaros 2005; Ryde & Svensson 2002; Ryde & Pe'er 
2009). The thermal component is expected to form a nar- 
row energy peak in the sub-MeV range, while a broad 
non-thermal spectrum describes the wide-band emission 
of GRBs. In Ryde & Pe'er (2009), the observed spectral 
time evolution of prompt emissions is described with a 
decrease in temperature of the thermal component in this 
thermal regime. 

As previous papers have argued, spectral evolution pro- 
vides a new perspective for understanding the physics 
in emission regions. The time scales of brightening and 
dimming are determined by several factors, such as the 
size and speed of the emission region and the acceler- 
ation/deceleration mechanism of electrons. The longest 
time scale among them dominates the observed light curve 
in the rise or decay phase, and its characteristic spectral 
evolution is observed in each phase. In case that the spec- 
trum exhibits a simple power-law shape, no spectral evo- 
lution is expected by variation of the emission region size 
or beaming effect alone. Even with the "curvature effect" , 
which accounts for the delayed photon emission from the 
high-latitude portion of the shell emission, it requires only 
a self-similar temporal decay, such as F v oc t~ 2 ~ a v~ a 
(Fenimore et al. 1996; Kumar & Panaitescu 2000; Demer 
2004; Qin 2008). However, once a spectral break crosses 
the observation energy range by deceleration of the emis- 
sion region, it is possible to create an apparent spectral 
evolution. In this regard, it has been shown in Zhang et 
al., (2009) that one can reproduce an exponentially de- 
caying light curve with an apparent spectral evolution by 
using a model invoking a power-law spectrum with expo- 
nential cutoff. On the other hand, heating/acceleration or 
cooling/ deceleration of electrons naturally produce each of 
the characteristic spectral evolutions. The time scales are 
subject to the emission mechanism and the physical state 
in the emission region, such as temperature variation, and 
the time scale of shock acceleration, or radiation electron 
cooling. 

To reduce the number of parameters affecting the 
time evolution, here we focus on fast-rise-and-exponential- 
decay (FRED) like light curves, which are the most 



promising for investigating the radiation process sepa- 
rately from the geometrical effects of the emission regions. 
The relatively fast rise implies that the time scale of ge- 
ometrical variation is sufficiently short not to dominate 
the longer decaying process. The exponential decay time 
scale is thus expected to reflect the state evolution of the 
emission region, including the radiation cooling process. 
In this paper, we investigate luminosity-spectrum evolu- 
tion by using data collected with the wide-band all-sky 
monitor (WAM; Yamaoka et al. 2009) onboard Suzaku 
(Mitsuda et al. 2007), and carried out spectral evolution 
study using energy-resolved light curves and time-resolved 
spectra. 

2. WAM observations and data selection 

The Suzaku WAM is equipped with the largest effective 
area among all GRB detectors currently in orbit. The de- 
tection part of WAM consists of 20 BGO scintillator units 
surrounding four sides of the main well-type phoswitch 
counter units of the narrow field Hard X-ray Detector as 
their anti-coincidence shield counters. The WAM moni- 
tors 27T radian of sky in energies between 50 and 5000 keV. 
The geometrical area of 800 cm 2 per one side and thick 
high-Z material (Z e ff — 71 for BGO) realized the effec- 
tive area of 400 cm 2 even at the incident photon energy of 
1 MeV. The WAM triggers a GRB onboard by using the 
rising edge of the count rate by 5.5 rms of the previous 8 s 
background count rate, and stores 4 energy-resolved light 
curves for the period between 8 s (16 s before 20 March 
2006) before and 56 s (112 s) after the trigger. The WAM 
data are stored in two different formats, namely: "TH" 
data, which contains the 4 energy-resolved light curves 
with high time resolution of 1/64 (1/32) s covering 64 ( 
128) s, and "PH" data. The energy band of each TH 
light curve is labeled as TH0, TH1, TH2, or TH3, which 
nominally cover the 50-110 keV, 110-220 keV, 220-540 
keV, and 540-5000 keV bands, respectively. "PH" data 
with 54 energy bin spectra and a 1-s time resolution cov- 
ering the energy band between 50 keV and 5 MeV, which 
are continuously accumulated to provide background in- 
formation for the Hard X-ray Detector. Since the WAM 
has no imaging capabilities, we can determine the angle 
of incidence only for a GRB whose position is determined 
from simultaneous observation results from the interplan- 
etary network (IPN) and/or Swift/BAT or Fermi/GBM. 
By using the position (and subsequently the incident an- 
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Fig. 1. TH0 light curves observed from the GRBs listed in table 1. The light curves are evaluated with the FRED model described 
in § 3 



gle), we generate a precise energy response for the spectral 
analysis with the PH data. 

We selected 6 GRBs with 7 isolated FRED-shaped 
pulses according to the following criteria in order to in- 
vestigate the spectral evolution from all GRB events trig- 
gered onboard between August 2005 and December 2010. 
(1) Bright: events whose TH2 count rate as observed by 
the most-irradiated side of WAM is higher than 1000 c s _1 
at peak, (2) Long decay: events exhibiting 3-s time bins 
containing more than 100 c s" 1 in energy bands up to 
TH2, (3) Fast-rise: events exhibiting fast rise in less than 
1/3 time interval of that of decay, in order to discriminate 
possible geometrical effect from the decay phase, and (4) 
Isolated: events containing well-separated (without over- 
lapping) pcak(s) that are suitable for evaluating spectral 
evolutions from these light curves. 

In table 1, information on selected events is sum- 
marized: GRB ID, trigger time, position (if deter- 
mined), duration including 90% of all observed pho- 



tons (T90), and a reference to a simultaneous detec- 
tion or a reported Gamma-ray burst Coordinate Network 
Circular (GCNC) number. The T90 values arc taken 
from the table published by the WAM team 1 . Since 
GRB 081224 was observed with Fermi-GBM (Wilson- 
Hodge & Connaughton 2008) and with Swift-XRT (Page 
& Evans 2008) in the afterglow phase, the position 
was well determined. GRB 100707A was observed with 
Knous-Wind (Golcnctskii et al. 2010), Fermi-GBM/LAT 
(Wilson-Hodge & Foley 2010;Pclassa & Pcsce-Rollins 
2010) and the position was sufficiently determined for the 
purpose of this study with the IPN (Hurley et al. 2010). 
Only for the both two events, we were able to generate 
energy responses using the positions as determined. A 
spectral analysis for each of these objects is presented in 
section 4. We used XSPEC version 12.5.1 for spectral 
analysis (Arnoud et al. 1996). All quoted errors are the 

1 http:/ /www. astro. isas.ac.jp/suzaku/research/HXD- 
WAM /WAM-GRB / 
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Table 2. Decay power-law index for each TH band 



GRB ID 


fa 


Pi 


fa 




(x 2 /d.o.f.) 


(xVd-o.f.) 


(x 2 /d.oi.) 


050924 


-0.82 ±0.05 


-0.91 ±0.05 


-0.86±0.10 




(71.5/16) 


(58.9/16) 


(40.5/12) 


060915(1) 


-0.87±0.15 


-1.38±0.18 


-0.88 ±0.24 




(34.7/8) 


(21.9/6) 


(3.1/6) 


060915(2) 


-1.06 ±0.07 


-1.03 ±0.18 


-1.30 ±0.63 




(26.8/18) 


(4.14/6) 


(0.14/2) 


060922 


-0.94 ±0.05 


-1.22 ± 0.05 


-1.39±0.15 




(66.0/16) 


(93.2/14) 


(9.80/6) 


061101 


-1.49 ±0.04 


-1.38±0.11 


-1.85 ±0.62 




(171/12) 


(13.6/6) 


(8.66/2) 


081224 


-0.87 ±0.06 


-1.06 ±0.07 


-1.12±0.15 




(30.1/14) 


(36.6/11) 


(8.99/6) 


100707A 


-0.84 ±005 


-0.97 ±0.08 


-1.33 ±0.27 




(45.6/9) 


(3.01/6) 


(9.3/6) 



The derived time constant f3 n (f(t) oc t°) is shown for each THn band 
(sec text). 



90% confidence level. 

3. Light curves 

In this section, we examine spectral evolution of the se- 
lected 7 FRED-shaped pulses using energy-resolved light 
curve (TH data). Figure 1 shows all the TH0 light curves 
from the selected events, including the background. Here 
we employed light curves binned to 0.5-s resolution in or- 
der to evaluate possible trend of spectral variability with 
sufficient statistics. We evaluated the averaged back- 
ground count rate (C) by fitting each light curve (f(t)) 
in TH0 - 3 separately, as shown in figure 1 for TH0, with 
the following function of t: 

{ C (t<t s ) 
f(t)=\ i^t Po + C (**<*<*p) 
[ Poe-^+C (ip<i), 

where t s and t p are the start time of the rising and the 
peak time; further, P and r are the peak count rate 
and the decay time constant, respectively. Since the light 
curves of GRB 061101 exhibit some activity before the 
FRED like peak, we omitted the time region before the 
peak from the fitting procedure. These evaluation with all 
the parameters left free shows that most light curves rise 
within 1 or 2 s but do not provide satisfactory statistics 
to discuss energy dependence. 

Then, we concentrated on the decaying portions after 
we subtracted the background count rate evaluated above. 
We omitted both time bin including the t p evaluated for 
TH0, which exhibits the latest time in each pulse, and 
the next time bin to avoid the possible effect by the rise 
portion. The end of the time regions were limited by con- 
sidering statistical and systematical uncertainty. Since 
the typical rms of the evaluated background countrate is 
~ 33 count s _1 , we omitted the time regions the lightcurve 



below 100 count s _1 to avoid possible background fluctu- 
ation. Consequently, all TH3 light curves were omitted 
from the following analysis. 

We firstly examined a power-law function fit for these 
background subtracted decay light curves. We set the 
time origin at t p in this power-law fit. However, none of 
the selected light curves were represented by a power-law 
function, as we see in table 2 summarizing the best fit de- 
cay index and the chi-squared value. Then we employed 
exponential decay functions to evaluate decaying time 
constant for each TH lightcurve. We left the paramter 
Pq free here. The exponential decay model succeeded in 
fitting all of the energy-divided light curves as shown in 
figure 1. Note we carefully ignored small (~ 200 c s _1 
in each band) post FRED activity in GRB 100707A to 
obtain successful fit in the estimation of time constant. 
Table 3 summarizes the best fit values and chi-squared 
values. 

Interestingly, as we see in table 3 and figure 2, all events 
exhibit shorter time scales in the higher energy bands. 
In order to evaluate this trend quantitatively, we fitted 
the derived time constants of decays with a power-law 
function (r = NE 1 ) and summarized the obtained best- 
fit value of 7 within the 90 % confidence level in table 3. 
The averaged value and its 90 % confidence level of 7 is 
-(0.34 ±0.12) (x 2 /d.o.f.= 3.7/6). It is noteworthy that 
the values of 7 as determined for 5 of the 7 FRED peaks 
are consistent with the average value of —1/2, which is 
the value expected for synchrotron or inverse- Compton 
cooling, while the rest 2 FRED peaks rejected. 

4. Time-resolved spectrum 

As mentioned in section 2, the positions of GRB 081224 
and GRB 100707A were determined with an accuracy suf- 
ficient to generate reliable energy response matrices. In 
this section, we analyze the respective time-resolved spec- 
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Fig. 2. FRED exponential decay features of the selected GRBs. The background subtracted light curves binned to 0.5-s resolution 
for THO, TH1, and TH2 bands are shown in each top, middle and bottom panels. The best-fit results are shown as solid lines and 
the end of evaluated time regions are indicated in dashed lines. 
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Table 3. Exponential decay time constants for each TH band and the energy index 









T 2 (S) 


7 1 




IX /a.o.i.j 


(v 2 /Hnfl 


(v 2 /dnfl 
IX / U.O.I. J 


(v 2 i (r\ n f — 11 

IX J (u-O.i. — l) 


At AHO A 

050924 


O AO 1 A 1 

3.93 ± 0.31 


O /I O 1 A O E 

3.43 ± 0.25 


O 1 n I A A A 

3.19 ± 0.44 


-0.14I 016 




fi q /i /i 


^ID.o/ ID J 


(lo. f / 1Z J 




060915(1) 


2.92 ± 0.53 


1.68 ± 0.25 


3.12 ± 1.01 


n a o+O 47 


(Zl.O/ 8J 


(14.7/6) 


(O.OZ 1 V) 


/q OOA 


060915(2) 


3.30 ±0.31 


2.39 ± 0.46 


1.46 ± 0.92 


^ r-r-0 44 

-0-45l^| 




(zo.4/ 10 J 


(i.oy/o j 


(U.47/2) 


(^U.UD ) 


060922 


3.58 ±0.31 


2.48 ±0.15 


1.60 ±0.21 






(27.1/16) 


(29.9/14) 


(5.05/6) 


(0.16) 


061101 


1.83 ±0.07 


1.67±0.15 


0.96 ±0.35 


-0 22+ 16 

u - zz -0.18 




(7.25/12) 


(3.20/6) 


(7.13/2) 


(1.60) 


081224 


3.56 ±0.32 


2.60 ±0.21 


1.93 ±0.31 


-0.40 ±0.18 




(10.9/14) 


(16.9/11) 


(2.93/6) 


(0.004) 


100707A 


2.97±0.20 


2.23 ±0.22 


1.49 ±0.40 


-0 40+ 19 




(1.41/9) 


(10.1/6) 


(16.0/6) 


(0.13) 



The derived time constant r n (f(t) oc e~ f / Tn ) is shown for each THn band (see text), 
t The energy index of power-law function (r oc E 7 ) to evaluate the time constant T n . 
Note that the degree of freedom is unity since the evaluation is carried out for the 
three energy divided time constants (see text). 





Fig. 3. The time-resolved folded spectra of GRB 081224 in the four second decay phases. Each 1-s integrated spectrum is fitted 
with a power law combined with a blackbody radiation model. The dotted lines represent the two components, while the solid 
histogram represents the sum of the two components. Ratio to the best fit models of BBPL and GRBM are shown in each middle 
and lower panels, respectively. 
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tra to examine the spectral behavior in detail. Exposure 
time of each spectrum is 1-s which is the time resolution 
of PH data of WAM. The time region here we present is 
limited to the time region that we showed in TH2 light 
curves in figure 2. In order to avoid uncertainty in gen- 
erating the energy response, we ignored the bins below 
150 keV in the following spectral analysis (Yamaoka et al. 
2009), while we used the upper most energy bin cover- 
ing up to 500 to 520 keV depending on the detector gain 
at each observation as provided by the WAM calibration 
team. 

Figure 3 shows four 1-s time-resolved spectra from 
GRB 081224 in the decay phase. We fit the four time- 
resolved spectra simultaneously to evaluate time trend 
of spectral parameters. Here we examined the spec- 
tra with a single power-law model (PL), an exponen- 
tial cut-off power-law model (CUT), GRBM (section 1) 
or a blackbody radiation with a powcr-law-shaped non- 
thermal emission model (BBPL). While PL failed to de- 
scribe the spectra with unacceptable chi-squared value of 
X 2 /d.o.f. = 527.9/59, CUT gave marginal fit with 68.6/55. 
On the other hand, the either of other two models - 
GRBM and BBPL provided significant improvement in 
fitting the data as are shown in figure 3. 

Fallowings are the notations we employed for the 
GRBM. 



/(*) = 




lOOkoV/ 
(ai— a 2 )Eo 
lOOkcV 
E \ Q 2 



(E< (a 1 -a 2 )E) ) 



. lOOkcV 

Figure 4 shows the time history of derived best-fit param- 
eters. We fitted the spectra in two steps. First we left 
all parameters free at the fitting and obtained the upper 
four panels shown in figure 4. Although these parame- 
ters describe the spectra well (x 2 / d.o.f. = 42.3/48), we 
cannot see any trends in the derived spectral parameters. 
Then, in order to evaluate the spectral evolution quanti- 
tatively avoiding the spectral parameter coupling, we tied 
the two photon indices of the four spectra. The obtained 
best fit values of lower and higher energy band photon 
indices in the simultaneous fitting were a\ = — (0.08lg 3 q) 
and a.2 = ~ (3-13lg'2g) with the chi-squared value of \ 2 1 
d.o.f.= 58.7/54. The other parameters, left free for the 
four spectra, normalization at 100 keV (No) and the char- 
acteristic break energy (Eq) are shown in bottom two 
panels of figure 4. The derived time trend of each of 
two parameters is well represented by a power-law func- 
tion of time, N oc t +1 - 26±0 - 30 ( x ^/d.ol .= 1.34/2) and 
E oc f-(°- 72 ±o.5i) (0.06/2). However, it is noticed that 
the trend of increasing of Nq and the hint of high energy 
photon excess in the ratio shown in figure 3 (d) . The for- 
mer suggests an artifact at describing "soft excess" with 
"steepening power-law with increasing normalization" or 
"lowering Eq with increasing normalization" as we see in 
figure 4. With the latter, suggesting high energy excess 
in the faintest phase, these features hint us a BBPL with 
fainting blackbody component describes the low and high 
energy excess better. 



TiiiiL' from trigger is) 



Time from trigger (s| 

Fig. 4. The time history of derived best-fit parameters from 
the time-resolved spectra of GRB 081224 as fitted with 
GRBM. The upper four panels show the results at four free 
parameter fit, while lower two panels show the derived best 
fit parameters of normalization and the characteristic break 
energy (i?o) f° r each second. In the latter case, the photon in- 
dices (011,012) below and above Eo are tied for all the spectra 
(see text). The normalization constant at 100 keV is in units 



of photons s 



^kcV~ 



and Eq is in keV, respectively. 



Time from [rigger (> I 

Fig. 5. The time history of derived best-fit parameters from 
the time-resolved spectra of GRB 081224 as fitted with a 
power law combined with a blackbody radiation model. From 
the top most to bottom panels show the time evolution of pho- 
ton index, 0.2 — lMeV flux of the power-law component in unit 
of erg s _1 cm" 2 , kT of the black body component in keV, and 
the normalization of the black body component in Lsi/D^q 
where L51 is the source luminosity in units of 10 51 erg s 
and Dio is the luminosity distance to the source in units of 
10 Gpc. The dashed line in each panel represents the best fit 
power-law function, whose best fit parameters are described 
in text. 
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Table 4. Time trends of best fit parameters for the time resolved spectra from GRB 081224 





BBPL model 






GR.BM 




parameter 


At 




parameter 


At 




a 

F(0.2 - lMeV) erg s" 
kT keV 

WD?o* 


o i 7 +().2t) 
Zli ' -0.28 

1 cm- 2 (1.3 ±0.3) x 10~ 6 
2101S 
70±1S 


(constant) 
-(0-43i°J 8 7 ) 

-(i.«)±g:9g) 


a\ 
a 2 
E keV 
N 


-(0-08^;-) 

-(3-i3±S:IS) 

OOU -260 

0.019 ±0.006 





-(0.72lg;«) 
±1.26 ±0.30 



t The time trend of each parameter are described in a power-law function At v . 

' Lsi/D^ where L51 is the source luminosity in units of 10 51 erg s — 1 and D10 is the luminosity distance to the source in units of 
10 Gpc. See text for the other paramters. 




200 500 1000 2000 5000 200 500 1000 2000 5000 

Energy (keV) Energy (keV) 




Fig. 6. The time-resolved folded spectra of GRB 100707A in each four second of the decay phase. Each 1-s integrated spectrum is 
fitted with a GRBM, which is shown in solid histograms. Ratios to the best fit model are shown in each lower panel. 



The BBPL succeeded to describe the time resolved spec- 
tra with x 2 / d.o.f. = 56.0/48. The parameters of the 
derived best fit with BBPL are shown in figure 5, which 
implies a relatively constant power-law component with 
a gradually cooling blackbody component. The photon 
indices of the derived best fit are consistent with a con- 
stant value of 2.19±g;i| (x 2 /d.o.f.= 1.16/4), while the 200- 
1000 keV flux shows an extremely small variation around 
the average value of (1.22 ±0.29) x 10~ 6 erg s _1 enr 2 
(x 2 /d.o.f.= 10.2/4). Interestingly, the derived behavior of 
kT is well reproduced by a power-law function of time with 
an index of -(0.43±g;^) (x 2 /d.o.f.= 0.02/2), whereas in 
the blackbody flux a monotonous power-law decay in the 
last 3 s with an index of -(0.47 ±0.33 (x 2 /d.o.f.= 8.3/2) 
although we see obvious deviation at the first bin. All 
derived time trend of the spectral parameters are summa- 
rized in table 4. Here, the results indicate that the energy 
dependence of the decay time constants in the energy- 
resolved light curves should be attributed to the cooling 
of the blackbody component qualitatively; however, from 
a quantitative point of view, the energy dependence is 



obscured by a fairly constant non-thermal component. 

The four 1-s time-resolved spectra from the decay phase 
in GRB 100707A exhibit no significant preference for the 
four and two parameter fit of GRBM (x 2 /d.o.f.= 62.4/48 
and = 62.8/54, respectively) and BBPL (x 2 /d.o.f.= 
60.6/48). However we cannot find any time trends with 
four nor two free parameters for each spectrum fit, be- 
cause of the large uncertainty of the best fit parameters. 
Moreover, similar chi-squared values for different d.o.f. 
suggests reducing free parameters. We examined single 
parameter fits by choosing the turn-over energy Eq for 
GRBM and tied the normalization, the first and second 
spectral index (ai,«2)- The time evolution of the best- fit 
value of Eq are summarized in figure 7. The parameter of 
Eq exhibits a power-law-shaped decrease in time with an 
index of — (1.11"q'|) (x 2 /d.o.f.= 0.03/2) as summarized in 
figure 7 and table 5. 
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Fig. 7. The time history of derived best-fit break energies 
(E ) in kcV for the time-resolved spectra of GRB 100707A. 
The other three spectral parameters (01,02, and the normal- 
ization constant) are tied. The dotted line indicates the best 
fit power-law function (see text). 

Table 5. Time trends of best fit parameters for the time re- 
solved spectra from GRB 100707A 





GRBM 




parameter 






ax 
a 2 
E keV 


-(0.83™) 
-(2-52^) 





-(1.1 ±0.7) 


N 


0.28 ±0.03 






function Af 1 . 

5. Discussion 

In order to investigate the cooling/deceleration mecha- 
nism, we selected 7 well-separated asymmetric FRED-likc 
peaks from 6 GRBs observed with WAM in accordance to 
the selection criteria described in section 2. As described 
in section 3, all decays are described well with expo- 
nential decay function while they clearly rejected power- 
law decay. We examined all of the energy-resolved light 
curves systematically and qualitatively with a linear-rise- 
exponential-decay function. Next, we found the obtained 
decay time constants, which are well described with a 
power-law function of the energy. All energy indices are 
negative, which implies shorter time constant at higher en- 
ergy band. Although as yet we cannot propose a clear in- 
terpretation for the light curves exhibiting 7 = 0.34 ±0.12 
in average, we note that the 5 FRED peaks showing 
7 = —1/2 are consistent with the expected results for 
non-thermal synchrotron and/or inverse- Compton cooling 
(Rybicki & Lightman 1979). 

On the other hand, as mentioned in section 1, no ex- 
ponential decay can be obtained with a simple power-law 
spectrum produced in a decelerated region. Even if we 
consider the curvature effect for a cut-off power-law spec- 
trum, it has been shown in Zhang et al., (2009) that the 
light curve exhibits a power-law type of decay in propor- 
tion to i~( 2+a ) around the cutoff energy, where a denotes 
the power-law index of the source spectrum. Therefore, we 
conclude that neither the observed exponential decay nor 
its energy dependence as observed in the 7 selected GRBs 
are expected in the geometrical interpretation. Therefore, 
at least for the cases selected under the abovementioned 
criteria, the decays are not dominated by geometrical ef- 
fects. 

Before concluding the emission mechanism or the time 



evolution of the physical parameters of the emission re- 
gion, we investigated the time-spectral evolution from 
the two GRBs located with the aid of simultaneous ob- 
servations. The time-resolved spectra obtained from 
GRB 081224, exhibiting energy index of decay time con- 
stant 7 = —0.40 ±0.18 in the light curves, are well repro- 
duced with BBPL as well as GRBM. However, it is not 
easy for us to draw a clear view of what makes both bright- 
ening with lowering characteristic break energy (Eq). 

The observed exponential decay in the energy-resolved 
light curves can also be attributed to the time evolution of 
the thermal component, however, the energy dependence 
observed in the light curve was under estimated by the rel- 
atively constant non-thermal component. It is also note- 
worthy that the derived temperature of the blackbody ra- 

— CO 43+ ' 27 ') 

diation component decreases in proportion to t 1 ' -o.as> 
(table 4). In the case of blackbody radiation cooling, the 
temperature variation at the source rest frame is governed 
by the following equation. Here, we assume conservation 
of the number of electrons in the emission region, and 
thus the product of the number density of electrons and 
the volume of the emission region nV is constant. 

dE 



dt dt 



T = 



n 

3aS 
nVk 



-1/3 



.-1/3 



where E, t, T, and S denote the total thermal energy, 
time, temperature, and emission area, respectively. We 
naturally regard the temperature decay (e.g T oc tr 1 ^ 3 ) 
observed from the spectral variation in GRB 081224 as a 
result of blackbody radiation cooling and/or a increase in 
the emission area with time (e.g. S oc t 2 ). 

On the other hand, the time-resolved spectra obtained 
from the other GRB 100707A are well reproduced with 
GRBM, although they do not accept a thermal model. 
The photon indices of the lower and higher energy bands 
{0-1,0.2) are consistent with those expected from the stan- 
dard synchrotron regime. The decay index of Eq is also 
consistent, though the uncertainty is large, with the dis- 
tribution reported in Peng et al. (2009), where FRED 
samples are examined using the data of the Burst and 
Transient Source Experiment on the Compton Gamma 
Ray Observatory ((decay index of Eq) ~ —0.73 or —0.76 
with a standard deviation of 0.22). The light curve in- 
dicating 7 = — (0.40ioag) suggests a non-thermal syn- 
chrotron or inverse-Compton process with a cooling time 
scale of r = E^ 1 / 2 (Rybicki & Lightman 1979). Adopting 
the synchrotron and/or inverse-Compton cooling regime, 
we estimated the magnetic field strength {B) as a function 
of the cooling frequency {u c ,obs) an d the cooling time scale 
(^c,obs) in the observer frame. By employing the source 
redshift (z), the Doppler factor (<5 = r/(l — /3 c cos0), where 
the Lorentz factor, velocity in the unit of light speed, and 
the angle between the jet and line of sight are denoted as 
r, j3 c , and 9, respectively), and the Compton parameter 
Y, B is expressed as a modification of Eq. (2) presented 
in Kumar & McMahon (2008), 
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where the physical constants electron charge, electron 
mass, speed of light, and Thomson cross-section are de- 
noted as q, ro e , c, and <tt, respectively. This equation can 
be rewritten as, 

-1/3 



B = l 



l + z 



1/3 



S 

To 3 

-1/3 



1 + Y 



20 



V lOOkeV ) 
l + z 



:,obs 



2 

-2/3 



-2/3 



1/3 



s 

10 3 



3s 

-1/3 



milligauss,or 



1 + Y 



10 



-2/3 



V10 14 Hz 



-V3 /t 



3s 



-2/3 



milligauss, 



in cases where synchrotron or inverse-Compton emissions 
are observed in the observed gamma ray band, respec- 
tively. Thus, a magnetic field with a respective strength 
of a few milligauss or a few tens milligauss can be expected 
for the two cases. These values are significantly smaller 
than those estimated in previous works. It is mainly 
because we employed relatively large cooling time scale 
events according to our selection criteria. This may sug- 
gest that the GRBs we studied had small beaming factors 
(8). Regarding the possible seed photons for the inverse- 
Compton process, it is necessary to consider not only syn- 
chrotron photons produced in the source (jet) rest frame, 
but also field photons outside the synchrotron electron 
field (e.g., Tashiro et al. (2009)). In particular, it is note- 
worthy that the energy density of cosmic microwave back- 
ground (CMB) photons, (ucmb), irradiating the GRB jet 
(ucmb = 4.1 x 10 _13 (1 + z) 4 T 3 , where z and F are the 
source redshift and the Lorentz factor of the jet, respec- 
tively) exceeds the energy density of the magnetic field. 
Assuming z = 1 and possible Lorentz factor e.g. T = 100 
or 1000, we then obtain «cmb = 6.5 x 10~ 6 erg cm 



or = 6.5 x 10~ 3 erg cm~ 3 , while the energy density of 

the magnetic field becomes ub 
-3 „_ „, a n v 1 n-5 



S 2 /(8tt) = 4.0 x 10~ 8 
erg cm~" or ub = 4.0 x 10~ 5 erg cm" 3 for synchrotron 
gamma rays or inverse-Compton gamma rays, respec- 
tively. Since the luminosity of synchrotron radiation and 
that of a CMB-boosted inverse-Compton process are pro- 
portional to ub and mcmb, respectively, this requires that 
the luminosity of a CMB-boosted Compton process be 
much greater than that of synchrotron radiation in the 
extremely relativistic jets assumed for the GRB. We need 
more precise information on the bulk Lorentz factor to 
conclude which the dominant seed photon source is. 
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